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ABSTRACT 

The afterglow emission that follows gamma-ray bursts (GRBs) contains valuable information 
about the circumburst medium and, therefore, about the GRB progenitor Theoretical studies 
of GRB blast waves, however, are often limited to simple density profiles for the external 
medium (mostly constant density and power-law /? * ones). We argue that a large fraction 
of long-duration GRBs should take place in massive stellar clusters where the circumburst 
medium is much more complicated. As a case study, we simulate the propagation of a GRB 
blast wave in a medium shaped by the collision of the winds of O and Wolf-Rayet stars, the 
typical distance of which is li ~ 0.1 - 1 pc. Assuming a spherical blast wave, the afterglow 
light curve shows a flattening followed by a shallow break on a timescale from hours up to a 
week after the burst, which is a result of the propagation of the blast wave through the shocked 
wind region. If the blast wave is collimated, the jet break may, in some cases, become very 
pronounced with the post break decline of the light curve as steep as r^. Inverse Compton 
scattering of ultra-violet photons from the nearby star off energetic electrons in the blast wave 
leads to a bright ~GeV afterglow flare that may be detectable by Fermi. 
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1 INTRODUCTION 

The prompt phase of GRBs is followed by a long-lived afterglow 
emission. The afterglow is believed to be powered by shocks driven 
by the relativistic outflow into the circumburst medium. Assuming 
that these external shocks inject ultrarelativistic electrons into the 
downstream shocked fluid, synchrotron emission from these parti- 
cles can account for th e basic properties of a large number of t he 
afterglow observations jSari et aljl998l ; lwiiers & Galamall 19991) . 

Afterglow modeling can provide important information about, 
among other things, the density profile of the circumburst mate- 
rial, thus constraining the nature of the progenitor. So far, however, 
typically only very simplistic density profiles have been consid- 
ered, i.e., the external medium is assumed either to have constant 
density or to follow a power-law function of spherical distance R. 
Observations are inconclusive about the density profile, i.e., de- 
pending on the burst and the analysis more than one or none of 
the simple density profiles may account for observations (see, e.g. , 
Panaitescu & Kumar 2002; Piro et al. 2005; Starling et al. 20oi; 



Curran et all2003l ; ISchulze et all201 ll) . 



At least some long-duration GRBs have been convincingly 
shown to be associated with the death of massive s t ars of Wolf- 
Rayet type (Galama etalj |l998l ; IStanek et alj |2003| ; iHiorth et all 
[2003; Mazzali et al. 2003j). Because Wolf-Rayet stars are charac- 
terized by strong winds, a R^^ wind-like profile is, at first sight, a 



* E-mail: Petar.Mimica@uv.es 



natural choice to describe the density distribution surrounding the 
progenitor. However, for high enough density of the medium that 
confines the stellar wind, the termination shock of the wind may 
take place s ufficiently close to the progenitor to affect the afterglow 
light curve jWiiersl2q0ll;IChevalier et alj2004lRamirez-Ruiz et all 
l2005l ; |Pe'er & Wiiersll2006l)" 

The actual density profile which decelerates the GRB -driven 
blast wave could be much more complicated. Massive stars rarely 
form in isolation; they preferentially reside in dense stellar clusters 
where tens or even hundreds Wolf-Rayet and O stars are crowded 
on sub-pc scale regions (e.g. iMassev & Hunted 1 19981) . About one 
third of the Galactic Wolf-R ayet stars are located in several very 
massive stellar clusters (e.g., lFigeill2004l) . It is reasonable to ex- 
pect a fair fraction of GRBs to take place in such dense stellar en- 
vironment. The interactions of the strong stellar winds in stellar 
clusters complicate the medium that surrounds the GRB progeni- 
tor and, therefore, the blast wave evolution. Furthermore, nearby O 
stars provide a strong UV photon field that is up-spattered by elec- 
trons accelerated at the forward shock and may potentially result in 
bright GeV afterglow flaring (jOLannios 2O0S). 

As a first step towards studying the afterglow appearance in 
more realistic density profiles, we focus on a blast wave propagat- 
ing through a medium shaped by colliding stellar winds. Hydrody- 
namic simulations are used to study the profile from the collision 
of stellar winds. Follow-up relativistic hydrodynamic simulations 
are performed in order to study the blast wave propagation through 
such density profile (Section |2ll. The fluid dynamics are coupled 
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Figure 1. Contours of the logarithm of the gas density (in arbitrary units) 
in the vicinity of the progenitor (P) and the O star. The contour values are 
decreasing in steps of 0.25 as one moves away from the stars. The shocked 
stellar winds are located in the region between two thick arcs, which are 
the termination shocks for the progenitor and O star winds. The three lines 
show the cuts at angles 5, 10 and 20 degrees from the line joining the centers 
of the two stars, respectively. 



to a radiative transfer code to calculate the resulting synchrotron 
and inverse Compton emission (Sections[3]and|4j. We discuss our 
results in Section|5] 



2 A BLAST WAVE IN COMPLEX DENSITY PROFILES 

With the number of massive stars N, ~ 100 (ranging from a few 
tens to a few hundreds) crowded on a ^ 1 pc scale of a typical 
massive stellar cluster (similar to, e.g., Westerlund 1 , Arches, Quin- 
tuplet, or Center in the Galaxy), the mean distance between O stars 
is fif ~ Rc/nI^^ ~ 0. 1 pc. The blast wave driven by a ORB, still rel- 
ativistic at these distances, is expected to encounter density bumps 
while propagating in the cluster. We simulate several density pro- 
files that may be expected in a massive cluster and then study the 
blast wave propagation through them. 



2.1 Colliding winds 

In the young stellar cluster under consideration, the gas density is 
shaped by the interactions of the stellar winds. Since the most mas- 
sive (O and Wolf-Rayet) stars have the strongest winds, they are 
going to dominate these interactions. The stars are surrounded by 
regions of their freely expanding winds that are followed, further 
out, by regions of shocked gas-result of wind- wind collisions. For 
a given line of sight from the explosion center to the observer, the 
density is likely to be shaped by the few massive stars that happen 
to lie close to it. For A', stars randomly distributed in the cluster, 
the closest one to the line of sight of the observer is located at an 
angle 9 ~ (4/N,yi^ ~ O.lN^l'^ rad (where the A = lO'A^ notation 
is adopted). Encountering a star within an angle 6 ~ 10° from the 
line of sight is, therefore, the norm in the massive clusters under 
consideration. As long as the outflow is ultrarelativistic with a bulk 
Lorentz factor F » 1, the observer of the GRB afterglow probes 
the blast-medium interactions that take place within a narrow cone 
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Figure 2. Density profiles along the cuts at 6 = 5. 10 and 20 degrees from 
the line joining the centers of the two stars on Fig.[T] The radius R is mea- 
sured from the center of the progenitor and is normalized to the distance 
between the stars. 



of angle ~ 1 /F with respect to the line of sight. In the narrow cone 
which the observer can see the closest stellar encounter is the dom- 
inant one in determining the relevant density profile. 

As a first approach, we limit ourselves to a single massive star 
(O) located at distance < 1 pc from the explosion center (P) and 
at a modest angle f < 0.5 rad with respect to the line that connects 
the explosion center to the observer (Fig.[TJ. We consider the den- 
sity profile resulting from the collision of the winds of two stars 
with mass loss rates Mp = lO'^Majr^' and Mo = lO^^Moyr"', 
respectively. The Mp is typical for a Wolf-Rayet star assumed to 
be the GRB progenitor, and Mq is expected for a typical O star 
(also referred to as the 'companion' star). The stellar winds are as- 
sumed to be cold and have the same velocity i'„ = 1000 km s"' . We 
have performed a 2-dimensional (2D) axisymmetric hydrodynam- 
ical simulation of such wind-wind interaction assuming adiabatic 
behaviour of the gafl We hav e used the high-resolution shock- 
capturing scheme MRGENESIS jMimica et al. 2007 , 2009bJ to run 
the simulation until a steady state is established. The numerical res- 
olution is 800 and 2400 cells in x and y directions, respectively. 

In Fig. [T] we show the density contours from such interac- 
tion. The regions of freely expanding winds and shocked winds 
are clearly seen. Because of the choice of the same velocity \'„ for 
the two winds, ther e is no contact disco ntinuity separating the two 
shocked winds (see lStevens et al The shocked wind region 

lies closer (and bends toward to) the O star because Mp > Mq. 

In Fig.|2] we plot the density profile moving radially outwards 
from the Wolf-Rayet star and for different angles with respect to the 
O star. The density profiles are averaged over a cone with opening 
angle of 5 degrees. One can see the R - density profile from the 
progenitor, a narrow region of enhanced density due to the shocked 
colliding winds and, at larger distance the freely expanding wind of 
the companion. These profiles are used as the density of the external 
medium through which the GRB blast wave propagates. 



' IStevens et alj ) 19921) and, more recently, Ivan Marie et alj )201 ih have 
studied the stellar wind interaction in detail and showed that for the dis- 
tances f/ > 0. 1 pc of interest here, radiative coohng in the shocked regions 
is negligible, thus justifying the assumption that the gas is adiabatic. 



GRB afterglows from realistic density profiles 3 



2.2 Blast wave propagation 

For characteristic distance d > 10" cm between the stars con- 
sidered here, the blast w ave has evolved into a self-similar 
( iBlandford & McKeelll976h stage before encountering the shocked 
wind region. As long as the bulk Lorentz factor of the shocked fluid 
is larger than the inverse opening angle of the jet: F^Sj' , 2D effects 
(e.g. lateral jet spreading) can be safely neglected. In this regime, 
the blast wave can be studied assuming a spherically symmetric 
evolution. We follow the evolution of the blast wave wit h ID rela- 
tivistic hydrodynamic simulations dMimica et a"il|2009lj ). We cau- 
tion, however, that by the end of our simulations, the fluid may 
decelerate to F ~ several and lateral spreading effects (not taken 
into account by our simulations) may be modestly important for 
9j ~ 0.1; a. value commonly inferred in GRBs (Frail et al. 2001). 
This point is discussed further in Section 3.1. 

Having specified the density profile for the external medium, 
we only need to choose the (isotropic equivalent) energy of the blast 
E for the self-similar initial conditions to be well defined. We con- 
sider a rather powerful GRB of E = 10^* erg. The distance between 
the stars is set lo d = 2.2 x 10'^ cm. The blast wave dynamics is 
followed with relativistic h ydrodynamic simulati ons using the code 
MRGENESIS described in lMimica etakl ( l2009bl) . The simulations 
have been performed in spherical symmetry with a numerical reso- 
lution of 16000 cells in the blast wave. We have generated five blast 
wave models: 

• three models with external medium density profiles shown in 
Fig. [2] N05 (cut along a line 5 degrees off the line joining centers 
of two stars), NIO (10 degrees) and A'20 (20 degrees) 

• a simulation in an external medium with a wind profile 
throughout: model WP 

• a simulation in an external medium with a wind profile out to 
the wind termination shock (TS) at 1.6 x 10"* cm (we assume a 
density jump of factor 4 at the TS and a constant density medium 
afterwards): model TS 

In Fig. [3] we show the Lorentz factor of the fluid at the forward 
shock as function of distance from the center of the explosion. The 
initial evolution of all models is that of a self-similar blast in a R^^ 
wind profil e where the Lorentz factor j ust behind the forward shock 
F oc i lBlandford & McKee|[l976l) . For models N05, NW and 
A'20 the self-similar evolution holds until distance R ~ d where the 
blast encounters the region of the shocked winds. At this stage, the 
enhancement of the density causes the forward shock to decelerate 
and a weak reverse shock launches into the blast. On exit from the 
shocked wind region, the density drops and the blast rarefies in the 
front. The density profile has a second bump due to the approach 
to the companion star. This results in one more drop in the Lorentz 
factor F, the depth of which depends on the angle 6. These non 
self-similar stages of the interaction can only be studied in detail 
using numerical simulations. At still larger distance the blast wave 
relaxes to a self-similar solution determined by the R^^ density pro- 
file ofthe_com£anmn_^to. to model WP the blast wave follows 
the IBlandford & McKe3 ( Il976h evolution, while in the model TS 
it decelerates much faster after encountering the constant density 
medium behind the termination shock (e.g. van Eerten et al. 2009). 



3 EMISSION 

After the hydro dynamical simulations have been performed we use 
the SPEV code l lMimica et al ] |2009ah to compute the light curves. 
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Figure 3. Lorentz factor of the shocked fluid behind the forward shock as 
a function of the distance from the progenitor for a blast wave propagating 
at an angle = 5, 10 and 20 degrees (thick full, thick dashed and thin full 
lines, respectively) from the line joining centers of the two stars (see also 
Fig.[T). Also sho wn is the evolution of the Lo rentz factor in a wind profile 
(dotted line), the IBlandford & McKed h976ll analytic solution for a wind 
profile (thick gray line closely following the dotted line) and the evolution 
in a wind profile with the termination shock (thin dashed line). The distance 
is normalized to the separation of the two stars. 



The details of how SPEV is applied to ca lculate the afterglow s 
from GRBs can be found in the section 3 o f iMimica et aP tOld) . 
We consider synchrotron emission from the shock-accelerated elec- 
trons (Sect. 3.1) and external inverse Compton scattering of the 
photon field from the companion star off the same electrons (Sect. 
3.2). In this paper we ignore the synchrotron self Compton pro- 
cess. Sect. 3 contains qualitative discussion and analytical estimates 
while the numerical results are presented in Sect. 4. 

3.1 Syndirotron emission 

During the initial self-similar stage of deceleration in the wind of 
the progenitor of density p = A/R', with A = Mp/4m'„ = 5 x 
10"A, gr-cm"' where A, = M_5V^'g, the bulk Lorentz factor of the 
fluid just behind the shock is 



9E 



l6nARc^ 



60eI''r:'''a-J'\ 



(1) 



The observer time evolves as t^bs = R/2T^c ^ 500R\jA,E^^ s (ne- 
glecting cosmological redshift effects). For the distance d ~ 3x lO" 
cm and the energy E ~ 10^^ erg, it takes typically a day for the blast 
wave to reach the shocked wind region (this timescale can range 
from hours to a week, depending on the various parameters). 

We make standard assumptions in calculating the s ynchrotron 
emissi on from shocked fluid in the blast wave following lSari et al] 
l ll998l) . i.e., we assume that a fraction e^, and £b of the dissipated 
energy goes into accelerating electrons into a power-law distri- 
bution with index p and amplifying magnetic field, respectively. 
Throughout this paper we fix = 0. 1, eg = 0.005, and p = 2.5. 
The light curve at the initial self-similar stage is that calculated by 
IChevalier&Lil ( l2000l) that focus on a blast propagating in a 
wind profile. At a time fobs ~ d/T^c the blast wave encounters the 
region of the shocked winds, which encounter causes a flattening 
in the light curve (Fig.|4j. This feature has already been discussed 
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in the case of a wind with a termination shock jNakar & GranotI 
I2OO7I : Ivan Eerten et alJlioO^ . In our setup, the blast wave crosses 
the shocked wind region on a short timescale and transits to a less 
dense wind of the companion. This transition, a distinct character- 
istic of colliding stellar winds, leads to a steeper decline of the light 
curve. More light curves and a discussion are presented in the next 
section. 

In our ID, spherically symmetric approximation, even the 
sharp features of the external medium (e.g. shocks, density jumps) 
result in smooth changes on the afterglow emission. This effect, 
result (at least in part) of the large lat eral extent of the emit - 
ting region, has been studi ed in detail by iNakar & Grano^ | |2007|) 
and Ivan Eerten et alj ( |2009|) . However, there is good evidence that 
GRB jets are co llimated with opening angles of Gj ~ 0.1 (e.g. 
lFrailetal.ll200ll) . Here we show that departures from spherical 
symmetry combined with structured external media introduce in- 
teresting novel features to the afterglow light curves. 

In addition to modulations of the light curve because of exter- 
nal density profile, deviations from spherical symmetry are also ex- 
pected to affect the afterglow appearance. A 'jet break' in the light 
curve is expected to occur whe n F ~ 67' f or a smooth density pro- 
file (Rhoads 1999). Zhang & MacFadverj | |2009|) have shown that 
the main effect of the T < Sj' transition is a steepening in the light 
curve because of th e 'missing surface' emitting towards the ob- 
server (see, however fwvgoda et al.l201 ll) . With our ID simulations 
we cannot treat the transition exactly. However, we can easily in- 
clude the dominant geometric effect that contributes to the jet break 
by considering the emission taking place only within an angle 6j 
with respect to the observer. Such examples are shown in Fig. |5] 
and discussed in Section 4. 1 . 

3.2 Compton scattering of the photons from the companion 

When the blast wave approaches the companion, it is exposed to the 
dense UV field of the O star. Electrons accelerated at the shock will 
Compton up-scatter the stellar light to (typically) GeV energies. 
Such encounter leads to a GeV flare hours to days after the GRB 
(the time of closer approach to the star) that is potentially detectable 
with Fermi, [oianniosl ( I2OO8I) analytically derived the y-ray fluence 
from a blast wave sweeping past the photon field of an O star for 
constant density ISM. Here we repeat the derivation for the wind 
profile, so that a comparison of the analytical estimates to more 
accurate numeri cal results can b e done (Section |4.2| (. In addition 
to the model of I Giannios! | |2008|) . we consider the contribution of 
the ambient UV photon field of the cluster to the external inverse 
Compton (EIC) scattering of blast wave electrons. 

We consider an O star of L, ~ lO^Lg ~ 10^"L395 erg-s"' 
located at a distance d from the progenitor, and at an angle 6 with 
respect to the line of sight. The emission of the star is assumed to 
peak at e, = lOe.j eV (typical for an O star). Electrons accelerated 
at the forward shock have a lower cut-off of their distribution at 
Tmin - (e,/3)F(m„/m,) - lOOOe^.-ifi^f rf7g"-A:"-, where Eq.Q is 
used in the last step. The energy of the scattered photons in the 
central engine frame is 



10^' r 



: 2F- 



J f,-i 54 



GeV. 



(2) 



The detailed numerical results (Section [4.2t identify the expression 
Cic as a good estimate for the peak of the L,, spectrum of the EIC 
component. Note that for e^^ ^ {nigC^f- le, ^ 25/e,j GeV, the scat- 
tering takes place in the Klein-Nishina regime for all electrons and 
the last expression is not applicable. 
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Figure 4. R-band (v = 4.3 X lO'* Hz) light curves for the five numerical 
models. Their light curves are indistinguishable until ^ 0.5 days, at which 
point the light curve of all the models which contain a shock (i.e., WOS, 
A'lO, ^20 and TS) flattens. The TS light curve soon reaches the asymp- 
totic behavior expected of the blas t wave propagating int o a homogeneous 
external medium (see e.g.. Fig. 2 of lvan Eerten et alj2009l) . The light curves 
of the remaining thi'ee models show a more complex behaviour due to the 
crossing over a second wind termination shock. The crossing point depends 
on the angle at which the blast wave is propagating, which is reflected in 
the different late-time light curves for WOS, A'lO and ^20 models. 



The fluence of the EIC component can be estimated by con- 
sidering the fraction of the total energy carried by electrons e^E that 
is radiated away because of EIC. The closest approach to the com- 
panion star is ~ 6d with the energy density of photons (in the rest 
frame of the blast wave) being (/ph - T^L,/4n(9d)^c. The IC cool- 
ing timescale of an electron with 7,^1,, is /cool = 2 x lO' /ymmUph s. 
The energy density of external photons peaks while the blast travels 
distance dd/T (in the comoving frame of the blast). The residence 
time in the intense radiation field is, therefore, /^s = r9/Tc. Com- 
bining all these, we obtain the fluence of the IC component: 



■q£ = 5x 10^ 



e^idj^A, 



erg. 



(3) 



The detailed numerical results (see Section 14^ show that expres- 
sion for overestimates the fluence of the EIC component by a 
factor of, typically, ~ a few. Klein-Nishina corrections contribute 
mostly to the discrepancy. 
The flare peaks at time 



'peak 



: d/2T-c ^ 0.5rf,8A.£54' days. 



(4) 



Notethatforrf ~ 3x10" cm. A, <1, a close encounter withabright 
star leads to extraction of a large fraction of the electron energy in 
the blast through EIC scattering. The emission from such encounter 
may peak several hours after the burst. 



4 LIGHT CURVES 

In this section we present the light curves produced by the five mod- 
els introduced in the section [2!2l We first discuss the optical. X-ray 
light curves result of synchrotron emission, and then turn our atten- 
tion to the y-ray emission result of EIC. 
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Figure 5. Same as Fig.|4] but for models where various jet opening angles 
have been assumed. Shown are the light curves of the N05 model assuming 
a spherical blast- wave (thick line), and assuming a jet half-opening angle of 
5° and 3° (thick dashed and dotted lines, respectively). For comparison we 
show the light curve in a wind environment for the spherical outflow and a 
jet half-opening angle of 3° (thin full and dashed lines, respectively), and 
two curves showing the decline proportional to r"* and (gray full and 
dashed lines, respectively). 



4.1 Optical/X-ray emission 

Fig.|4]shows the optical light curves. Before reaching the progenitor 
wind termination shock the light curves are indistinguishable from 
the one corresponding to a self-similar blast wave propagating in 
a R^^ profile (model WP). From that point on the models begin to 
diverge depending on whether the blast wave continues propagating 
into a constant-density environment (TS ) or whether it eventually 
encounters the wind of the companion star (A'05, NIO and ^20). 
The latter three models differ in their light curves because the blast 
wave crosses the wind interaction zone and a companion stellar 
wind at different angles. 

In model A'05 the blast wave propagates closest to the com- 
panion star and it encounters the densest companion wind. There- 
fore it is slowed down more than A' 10, which in turn decelerates 
more than A'20 (see also Fig.[3l(- This is seen in Fig. |4] after 1^3 
days, where A'05 is the brightest of the three models, followed by 
A' 10 and N20. 

Fig. [5] shows the effect of the finite opening angle of the jet. 
Since we are simulating a one-dimensional spherical blast wave, 
we model a jet with a half-opening angle Oj by assuming no con- 
tribution to emission from the fluid at angle with respect to the line 
of sight > 9j. As expected, while F > fy' the jet collimation is not 
affecting the light curve. For a blast wave propagating in a pure 
wind-like profile, a rather smooth break appears in the light curve 
when F ~ The break can be much sharper if the transition to 
F < 6t' takes place when the blast wave reaches the shocked wind 
region. The rapid decline of the Lorentz factor of the blast wave at 
the density jump is not compensated by the increase of the emitting 
surface visible to the observer and the flux drops much faster than 
expected from a jet break in a smooth external medium. Comparing 
with the thick full and dashed gray lines (which show the temporal 
decline proportional to r* and r^, respectively), we see that the 
combination of the jet break and the interaction of the blast wave 
with a shocked wind environment can produce steep declines. 

Fig. E] shows also demonstrates how the break happens at pro- 
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Figure 6. 1 keV X-ray light curve for the models A'05, A'lO and A'20 (thick 
full, dashed and dotted lines, respectively). Shown are also two cases for the 
model A'05 if we assume a jet half-opening angle of 5° and 3° (thin dashed 
and dotted lines, respectively). 



gressively earlier times as 0^ is decreases (thick full, dashed and 
dotted lines, respectively). At later times the light curve decline 
becomes less steep due to the acceleration of the blast wave as it 
leaves the shocked region and encounters a less dense companion 
stellar wind. 

Fig.|6]shows the X-ray light curve for the models A'05, A'lO 
and A'20, as well as for the model A'05 assuming a small jet opening 
angle, to simulate the effect of a jet break. As can be seen, the effect 
of the traversal of the shocked wind leaves qualitatively similar but 
less pronounced imprints on the X-rays in comparison to the optical 
light curve (see Fig.|4]l. The effect of the jet break (result of a finite 
jet angle) is, as expected, also an adiabatic one (e.g., not due to 
electron cooling), appearing in simultaneously in the optical and 
X-rays. 

4.2 y-ray emission 

The analytical estimates of Section IT2l provide the dependence of 
the photon energy and fluence (e,^ and F,^, respectively) of the EIC 
emission on the various parameters. The actual numerical values of 
Eqs.|2]and[3]are, however, meant more as rough order-of-magnitude 
estimates then the accurate predictions. In this section we use the 
numerical simulations to calculate the EIC emission and calibrate 
the analytics. 

Part of the uncertainty of the analytical estimates is connected 
to the fact that they ignore the blast wave hydrodynamics of the 
shocked wind regions (instead, the blast wave is assumed to prop- 
agate in a freely propagating stellar wind). The numerical simula- 
tions are more accurate since they include the effect of the collid- 
ing winds on the blast dynamics and the exact calculation of the 
EIC cooling during the blast encounter (including Klein-Nishina 
effects) for a power-law injected electron distribution. We compute 
the EIC emission assuming a monochromatic external point-source 
of radiation (good approximation for the black body stellar emis- 
sion). To compute the emissivity we use the method described in 
the section 2.2. 1 of lMimical ( |2004|) FI 

^ This enables us to use the same piecewise power-law representation of the 
non-thermal electron distribution that SPEV uses to compute synchrotron 
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Figure 7. Bolometric (4 MeV - 4 TeV) light curves from external inverse 
Compton scattering for the models A'05, A' 10 and A'20 (full, dashed and 
dotted lines, respectively). The thick lines show the light curves assuming 
that the only source of photons is the companion star The numbers in bold 
show the energy of the spectral peak in GeV for radiation observed 1,2, and 
10 days after the burst. The thin lines show the emission when the diffuse 
cluster emission is taken into account as a soft-photon source (see text for 
details). The numbers in italic denote the spectral peak in GeV for this case. 



Fig. |7] shows the EIC bolometric light curves and the energy 
of the spectral peak (in GeV) for the models A'05, N\Q and ^20 
(thick full, dashed and dotted lines, respectively). The EIC emis- 
sion peaks at time t^^ak ~ 2 days in good agreement to Eq.|4] (using 
the value dig = 2.2) and then gradually declines. The duration of 
the high energy flare is 6t ~ fpeak- As expected, the EIC is brighter 
for the closer encounters (smaller S) between the line of sight and 
the companion star. At maximum of the EIC luminosity the peak 
of the Ly spectrum is £ 2 GeV independently of 6 in close agree- 
ment to Eq. [2] The spectral peak moves to lower energies as func- 
tion of time. We have verified that The EIC spectrum has a cutoff at 
~ 25 GeV because of the Klein-Nishina suppression. The numeri- 
cally calculated fluence is found to be a factor of ~ 3 less that the 
estimate in Eq.[3] 

The total number of y-ray photons emitted in the 0. 1-300 GeV 
energy range (where the effective area of Fermi-LAr peaks) is Ny ^ 
8 X 10'', 2 X 10"", 6 X 10™ for = 5M0°,20°, respectively. For 
an effective area of LAT of A^ff ~ lO'' cm^, the closest encounter 
can be detected out to a proper distance of rfp ^ {AgffNyjAnY^^ ^ 
2.5x10^^ cm or out to a modest redshift of z ~ 0.3 for this example. 

In addition to the photon field of the companion, the blast 
wave also encounters the diffuse emission of the stellar cluster. The 
latter is dominated by that of the most massive stars in the cluster. 
For A', massive stars of luminosity L, each distributed isotropically 
within a cluster of radius R^, the ambient UV photon field in the 
cluster is U\}\ ~ 3N,L,/47tRIc. Fig.[7]shows the total EIC emission 
coming for the scattering of the photon field of the companion and 
the ambient cluster emission for A', = 30, L, = 10^'"' erg-s"', and 
Rc = 5 X lO'* cm (thin lines). Note that the ambient photon field 
is the dominant source of soft photons for the more distant encoun- 
ters (with 6 = 10, 20 deg) while it makes a modest contribution to 
the 8 = 5 deg example. The peak of the EIC component is broader 



emission, resulting in a considerable speed-up of the EIC emissivity calcu- 
lation. 



in this case. Including both the companion and the diffuse sources 
of soft photons, the total number of y-ray photons emitted in the 
0.1-300 GeV energy range increases to A'^ ^ 1.4 x 10'^, 8 x 10", 
6 X 10" for e = 5°, 10°, 20°, respectively. In this case the emission 
can be detected out to rf„ ^ 3.3 x 10"' cm, or out to z ~ 0.3 - 0.4. 



5 DISCUSSION 

If all, or at least the majority of, long duration GRBs come from the 
death of massive stars, a fair fraction of them should take place in- 
side young, massive stellar clusters. In such crowded environments 
the wind of the progenitor star terminates on sub-pc scales (mainly 
due to interactions with other stellar winds). The wind- wind inter- 
action regions are characterized by shocks, contact discontinuities 
and other sharp features in the density profile. The blast wave, result 
of the GRB jet interacting with the external medium, has a bumpy 
ride in such environments. In this work we search for the character- 
istic observational signatures of such blast wave evolution, focusing 
on the profile shaped by the collision of the wind of the progenitor 
with that of a nearby massive star or 'companion'. 



5.1 Optical and X-ray emission 

As long as the Lorentz factor of the blast F » 1 and F > fly', one 
can study the basic features of the blast-external medium interac- 
tion assuming spherical symmetry and considering different lines 
of sight to the observer. We find that, when the blast wave enters 
the shocked wind region, the optical light curve at first flattens, fol- 
lowed by a steeper decline at later times (Fig.|4}. When the charac- 
teristic synchrotron frequency crosses the observed band, the mod- 
ulations of the light curve are enhanced. Qualitatively similar (but 
much less pronounced effects) are in general expected in the X- 
rays as well (Fig.|6](. O ur results are in agreem ent with the generic 
arguments presented in lNakar & GranoJ l l2007h . who show that rel- 
ativistic blast waves in spherical symmetry do not show sharp fea- 
tures in their light curve even for step-function changes of the ex- 
ternal density as function of radius. 

Observations with dense sampling until late times and good 
signal-to-noise are required to detect the relatively smooth changes 
in the light curves. One recent examp le is a well studi ed opti- 
cal. X-ray afterglow of GRB 0804 13B l lFilgas et aklbouh . which 
shows a characteristic flattening in the optical followed by a 
steeper decline in both optical and X-ray wavelengths that ap- 
pears to be in agreement with the expectati ons from our model . 
Further examples include GRBs 21004 (see Lazzati et alj l l2002h 
and references there in), 080710 jKruhler et ai] 2009 j) , 080129 
( iGreiner et al]|2009ah and 080913 l lGreiner et aU 2009b). Finally, 
there is an interesting GRB 071031 jKruhler et al.,.2009h) . which 
shows several bumps in the optical, which might be due to more 
than one star influencing the afterglow (see, however, the discus- 
sion in Sec. |5.3| l. 

For a wide range of parameters, the blast wave is expected 
to enter the shocked wind region when F ~ (?t'. Strictly speak- 
ing, our ID calculations are not directly applicable to this regime 
since they ignore lateral spreading of the blast. Nevertheless, recent 
2D simulations performed by Zhang & MacFadven ( 2009) indicate 
that lateral spreading effects are small at this stage and that the 'jet 
breaks' are mainly a geometric effect. We show that the combina- 
tion of jet coUimation and structured external medium may lead to 
the rather sharp features in the light curve, such as a fast decline of 
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the light curve as steep as F oc f which is not expected for the 
evolution in a smooth external medium (Fig.|5]l- 

5.2 y-ray emission 

Another characteristic signature of the companion-blast interaction 
is an inverse-Compton powered, GeV flare (Giannios 2008). We 
show that, under optimistic conditions, a large fraction of the en- 
ergy deposited on the shocked electrons is radiated away during 
the encounter with the dense UV photon field of the massive star. 
Such a flare at the energies of tens of GeV may account e.g., for the 
emission of the GRB 940217 observed hours after the burst with 
EGRET ( lHurievlll994l) . The LAf instrument on-board to the Fermi 
satellite is capable of detecting such flaring for bursts out to mod- 
est redshift z < 1. Depending of the distance of the companion star, 
the ~ GeV emission will peak on timescales that range from a few 
hours up to a few weeks after the burst (Fig.|7]l. In some cases, the 
ambient UV photon field of the massive stars in the cluster may be 
the dominant source of the external inverse Compton component 
resulting in slower varying y-ray emission. A systematic search for 
delayed y-ray emission from the location of the burst on timescales 
of hours to weeks after the burst may be fruitful. The timing and 
photon-energy of such detection will provide invaluable informa- 
tion about the GRB environment. 

An additional promising source of seed photons to be inverse 
Compton scattered at the forward shock may be the infrared emis- 
sion from dust. Though its physical origin is unclear, dust is form- 
ing actively in th e stellar (or sho cked) wind region in many Wolf- 
Ray et stars (e.g., ICrowthej|2003l) . Up to ~ 10% of the stellar light 
can be reprocessed into the infrared that provides an isotropic bath 
of photons through which the blast propagates. Mid-to-near infra 
red photons can be scattered into multi TeV energies by the high- 
energy tail of the electron distribution while still in the Thomson 
scattering regime. Such emission, with a v ■ L,, peak at ~ TeV ener- 
gies, may be detectable by atmospheric Cherenkov telescopes pro- 
vided that the burst takes place at redshift z 0.5 (for the TeV pho- 
tons to avoid attenuation while propagating through the extragalac- 
tic background light). 

5.3 Other effects 

In this paper, we focused on the effect of a single massive star on 
the circumburst medium density in the cluster where the GRB takes 
place. This star is most relevant in the early afterglow stages since 
it lies the closest to the line of sight to the observer. Clearly, the 
collective effect of many stellar winds will dominate the density 
profile at larger distance (more than a few times the distance to 
the 'companion'). At these scales the external medium is expected 
to be very inhomogeneous, although on scales much smaller than 
those relevant for the blast wave emitting toward the observer. At 
late times, therefore, the afterglow light curve may resemble that 
expected from a constant density medium. The confining effect of 
the environment in the cluster can thus naturally limit the extend of 
the wind of the progenitor and lead to a constant-like density as is 
co mmonly suggested b y modeling of afterglow observations (e.g., 
see lSchulze et al.ll201lh . 
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